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ABSTRACT 

Galactic Center (GC) molecular cloud Sgr B2 is the best manifestation of an X-ray reflection nebula 
(XRN) reprocessing a past giant outburst from the supermassive black hole Sgr A*. Alternatively, 

Sgr B2 could be illuminated by low-energy cosmic ray electrons (LECRe) or protons (LECRp). In 
2013, NuSTAR for the first time resolved Sgr B2 hard X-ray emission on sub-arcminute scales. Two 
prominent features are detected above 10 keV - a newly emerging cloud GO.66—0.13 and the central 
90" radius region containing two compact cores Sgr B2(M) and Sgr B2(N) surrounded by diffuse 
emission. It is inconclusive whether the remaining level of Sgr B2 emission is still decreasing or 
has reached a constant background level. A decreasing Fe Ka emission can be best explained by 
XRN while a constant background emission can be best explained by LECRp. In the XRN scenario, 
the 3-79 keV Sgr B2 spectrum can well constrain the past Sgr A* outburst, resulting in an outburst 
spectrum with a peak luminosity of L3-79 keV ~ 5 x 10^® erg s“® derived from the maximum Compton- 
scattered continuum and the Fe Ka emission consistently. The XRN scenario is preferred by the 
fast variability of GO.66—0.13, which could be a molecular clump located in the Sgr B2 envelope 
reflecting the same Sgr A* outburst. In the LECRp scenario, we derived the required CR ion power 
dW/dt = (1 — 4) X 10®® erg s“® and the CR ionization rate Clt = (6 ~ 10) x 10“®® H“® s“®. The 
Sgr B2 background level X-ray emission will be a powerful tool to constrain GC CR population. 

Subject headings: Galaxyrcenter — X-rays: individual (Sgr B2) — X-rays: ISM — Molecular Clouds 


1. INTRODUCTION 


The Galactic center (GC) supermassive black hole 
Sagittarius A* (Sgr A*) is the closest such object and 
is thus an ideal target for in vestigation of galactic nu- 


clei and their activity cycles (Morris et al. 
Sgr A 


1999 Ponti 
is an underluminous black hole 


et al.||2013| ). 

with a bolometric luminosity about 10“® times th e Ed¬ 


dington luminosity for a 4 x lO ^Affn black hole (Chez 


et ai.||20()8 Gillessen et al.|[2009 ). Its current X-ray qui- 
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escent state, with a lu minosity of Lx ~ 10®® erg s 


-1 


a few times 10'®° erg 

s ® (e.g. 

Baganoff et al. 

2001 

Porquet et al. 

2008 

Nowak et a 

. 

21)12), during which 


lower than its Eddington luminosity. Hard X-rays up 
to ^ 79 ke V have also been det ected from the flares by 


NuSTAR (Barriere et al. 
experienced more substan 


2014). Whether it has ever 
lal increases in activity as ob¬ 


served in low-luminosity Active Galactic Nuclei (AGN) 
is still under discussion. 

Indication of such past activity of Sgr A* has come 
from Galactic center molecular clouds (GC MCs). In 


the populous Central Molecular Zone (CMZ, Morris & 


Serabyn| 1996), Sgr B2 is the densest and most mas- 
sive molecular cloud. It has complicated substructures 
including compact star-fo rming cores like Sgr B2(M), 
Sgr B2(N) and Sgr B2(S) (Etxaluze et al. 2013 and ref¬ 
erences therein). In an extremely simplihed picture, its 
density distribution can be described as a dense core with 
a radius of 0.15-0.3 pc (2"-4" assumi ng the Sgr B2 dis¬ 
tance of 7.9±0.8 kpc, Reid et al.|2009 ) and an H2 density 
of (3-9) X10® cm“®, surrounded by an envelope extended 
to ^ 5 pc (^ 2.2') with a density of 10'®-10® cm“® and 
a larger diffuse component reaching ^ 22 .5 pc with a 
roughly constant density of 10® cm“® (|Lis & Gold- 
smith||1990t |de Vicente et al.||1997|). 


Sgr B2 is the first CCAIC from which a strong 6.4 keV 
Fe Kq; line was discovered ( Koyama et al.|[l9M ). ASCA 
detected this significant line feature with an equivalent 
width of about 1 keV, which was later confirmed by 
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Chandra, Suzaku and XMM-Newton observations dM^ 


rakami et al.| 2001 Koyama et al.||2007 'i'errier et al. 
2010p. Time variability of the Sgr B2 Fe Ka line was 
revealed by years of monitoring by different instruments. 
The 6.4 keV line flux began declinin g in 2001 and de - 
creased by a factor of ~ 0.4 by 2005 (Inui et al.||2009|), 
and further decreased by a factor of 0.39±0.04 from 2005 


to 2009 ( Nobukawa et al.]|2011 ). The flux and morphol¬ 


ogy chanige of the X-ray emission is also observed from 
other GC molecular clouds. Interestingly, the clouds ex¬ 
hibit Fe Ka line variability in different ways, some rising, 
others are decreasing and an emission peak is detected 


emission fro m several GC molecular clouds (Yusef-Zadeh 
et al. 20071, it m eets challenges for most of the k'e Ka 
bright structures (Dogiel et al.||2013 20141. In the case 
of Sgr B2, the cosrnic ray electron energy required to pro¬ 
duce the cloud X-ray emission is as high as the bolometri c 
luminosity of the entire cloud (Revnivtsev et al. 2004). 
The derived metallicity of ZpZiQ ~ 3.1 is also much 
higher than current measurements of the GC metallic¬ 
ity, wh ich r anges from slightly h igher than solar (^un^ 


et al. 2008 Davies et al. 


2009) to twice solar (Giveon 


et al.|2002 ). Howev er, even it LFCRe is not the dominant 
process m the GC (Dogiel et al.| 2013) and particularly 

ini 


two years ( 

Muno et al. 

2007 

Ponti et al. 

2010 

Capelli 

et al.|2012 

Clavel et al. 

2UT3). 


rial is a natural explanation of the observed X-ray emis¬ 
sion. In the X-ray reflection nebula (XRN) model, K- 
shell photo-ionization and the subsequent fluorescence 
produce a strong Fe Ka line with an equivalent width 
EW> 1 keV, while a competition between Compton scat¬ 
tering of high energy photons and photoelectric absorp¬ 
tion of low energy photons gives ri se to a Compton re¬ 


in Sgr B2, we cannot exclude that, in specific regions, 
the LECRe process contributes to the background level 
of the Fe Ka emission. 

Though the GCMCs have been studied extensively be¬ 
low 10 keV, the investigation of the continuum emission 
extending beyond 10 keV has been limited. The hard 
X-rays from the S gr B2 region were f irst detected by 
GRANAT/kBU-Y’ ( Sunyaev et al.|1993|) and then by I/V- 


TEGRAL/lYASi in 2004 (iRevnivtsev et al. 


20041. In the 

XRN scenario, I'i'errier et 'al.|(|2U10l) derived an illuminat- 


flection hump around 20 — 

30 keV (Sunyaev et al. 111993 

Sunyaev & Churazov|1998| 

Koyama et ai.|1996|. An em- 


be an illuminating source, since no transient source has 
been sufficiently bright (L > 10^^ erg s“^) f or several 


ing source spectral index of I ~ 2 with combined XMM- 
Newton and INTEGRAL/IBIS spectra. Years of moni¬ 
toring with INTEGRAL/IBIS reveal that the hard X-ray 
emission decreases b y a factor of 0.4 from 2003 to 2009 
( Terrier et al.|^010| . Nevertheless, INTEGRAL/IBIS 


years since 1993 within or close to Sgr B2 (e.g. Revnivt- 
sev et al.|2004 Terrier et al.|2010| . Sgr A* was proposed 


to be the likely source illuminating Sgr B2. A major out¬ 
burst from Sgr A*, with a luminosity of a few 10^® erg s“^ 
lasting more than 10 years and ending a few hundred 
years 


was not able to resolye the hard X-ray emission. The 
accuracy of the hard X-ray luminosity was limited by 
the unknown distribution of the emission. 

Hard X-ray obseryations are crucial to better constrain 
the origin of the GCMC non-thermal X-ray emission. 
With unprecedented spatial and spectral resolution in 
the 10-79 keV band, NuSTAR obseryed GCMCs in 2013 


et al. 

1996 

Terrier et al. 

2010 

). 'I'his hypothesis is rein- as part ot the Galactic plane survey campaign, including 

lorcec 

by the discovery ol 

a superluminal h'e K echo trom the Sgr B2 region. I'br the hrst time, ISIuS'I'AK resolved 


the “Bridge”, which po ints to propagation of an eyent far 


away from the clouds (Ponti et al. 2010). Howeyer, the 
story gets complicated by ditterent Fe Ka line yariabil- 
ity detected from GCMCs, which cannot be explained 
with a single outburst from Sgr A*. The strong and fast 
yariation of the Fe Ka line flux in the “Bridge” region 
requires a two-year peaked outburst with luminosity of 
at least 10^® erg s“^, while the slower Fe Ka line behav¬ 
ior detected in other clouds suggests a second flare with 


a longer duration (Clavel et al. 2013). 

The propagatiori of cosmic-ray (GR) particles within 
the molecular clouds is an alternative explanation. 
Low-energy cosmic ray electrons (LECRe) and protons 
(LE CRp) can both produce hard X-rays and Fe Ka emis¬ 
sion ( Valinia et al.|2000 Dogiel et al.|2009 ). The Sgr B2 


Fe Ka hux variation time scale ot ~ 10 years is too short 
compared to the Coulomb cooling time of ~ 100 MeV 
protons invoked in the LECRp scenario, thus rules it out 
as a major cont ributor to the fast ch anging Fe Ka emis¬ 


sion of Sgr B2 (Terrier et al. 2010). However, LECRp 
could be a major contributor to the constant background 
level of the Sgr B2 Ee Ka emission, which could be 
detect able once the refle c ted X-ray em i ssion completely 
fades (iDogiel et al. 120091) . iDogiel et al.l (|2009|) estimated 


that LECRp could contribute to about 15% of the ob¬ 
served maximum Fe Ka flux obtained around 2000. In 
the LECRe scenario, a fast variation can be reproduced. 
Even if LECRe seems to successfully explain the X-ray 


its hard X-ray morphology and obtained a broadband 
spectrum from a single instrument. In Section 2, we in¬ 
troduce the observations and data reduction for NuSTAR 
and XMM-Newton data used for the following analysis. 
We present the morphology of the central 90" radius re¬ 
gion of Sgr B2 and the newly discovered cloud feature 
GO.66—0.13 in Section 3, their time variability in Sec¬ 
tion 4 and their spectroscopy in Section 5. Based on 
these observational results, we discuss their implication 
on both XRN and LECR scenarios in Section 6. 

2. OBSERVATION AND DATA REDUCTION 
2.1. NuSTAR Data 

The NuSTAR observatory oper ates in the broad X - 
ray energy band from 3 to 79 keV (Harrison et al.|2013). 
Sgr B2 was observed by NuSTAR in October 2013 in two 
25% overlapping pointings, with a total exposure time of 
293.7 ks (See Table 0. 

In both observations, the Sgr B2 region was imaged 
with the two co-aligned X-ray telescopes, with corre¬ 
sponding focal plane modules FPMA and FPMB, pro¬ 
viding an angular resolution of 58" Half Power Diame¬ 
ter (HPD) and 18" Eull Width Half Maximum (EWHM) 
over the 3-79 keV X-ray band, with a characteristic spec¬ 
tral resolution of 400 eV (EWHM) at 10 keV. The nom¬ 
inal reconstructed NuSTA R astrometry is accu rate to 
8" (90% confidence level, Harrison et al. 2013). The 
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data were reduced and analyzed using the NuSTAR Data 
Analysis Software NuSTARDAS v. 1.3.1. and HEA- 
SOFT V. 6.13, then hltered for periods of high instrumen¬ 
tal background due to South Atlantic Anomaly (SAA) 
passages and known bad/noisy detector pixels. 

The NuSTAR detectors are not completely shielded 
from incident X-rays that do not go through the op¬ 
tics, which is referred to as stray light. Bright X-ray 
sources within ~ 1-5 degrees of the NuSTAR field of 
view can significantly contaminate one or both of the de¬ 
tector planes. The contaminated detector pixels can be 
removed based on a numerical mo del that fully takes int o 
account the telescope geometry (Krivonos et al. 2014). 
Thus, to make a stray-light-free mosaic, we used this 
model to remove stray-light patterns from both FPMA 
and FPMB detectors by flagging the contaminated de¬ 
tector pixels as bad, when processing with the NuSTAR¬ 
DAS pipeline. As a result, stray-light from the the X- 
ray sources SLX 1744—299 and IE 1740.7—2942 were 
removed from FPMA detectors, and that from GX 3-1-1 
and SLX 1735—269 were removed from FPMB detectors. 
We then registered these stray-light-free images with the 
brightest point sources available in individual observa¬ 
tions. As a result, the astrometry is improved to ^ 4". 
Lastly we combined the exposure-corrected images in dif¬ 
ferent energy bands (Figure 1). 

On the other hand, removing the stray-light patterns 
is not necessary for spectral extraction. The stray-light 
background does not significantly change from one sky 
pointing to another, pr ovided that they are s eparated by 
no more than 10'-20' (Krivonos et al. 2014). Therefore, 


Table 1 

NuSTAR and XMM-Newton observations of Sgr B2. 


for the Sgr B2 core region detected in one observation, 
we extracted the background from the same detector re¬ 
gion in the other observation. This was our motivation 
for using two 25% overlapping sky pointings. This back¬ 
ground subtraction method has been applied to many 
sourc es detected by NuSTAR suffering from st ray-light 


(e.g. Zhang et al. 


2014 


Krivonos et al. 2014), and is 
proven to subtract mild stray-light contamination prop- 
erly. 

Due to the extremely bright stray-light contamination 
within 90" of Sgr B2 on FPMB, we used FPMA data only 
for the Sgr B2 core spectrum. From the selected data 
sets, we extracted source spectra from a circular region 
of 90" radius centered on Sgr B2 (R.A.=17'^47™20.4^ 
Decl.=—28^23™07.0®, J2000). In this way, we had two 
source-background pairs with the Sgr B2 core detected 
on FPMA in observation 40012018002 and 40012019001. 
We combined the two FPMA source spectra and their 
associated response files and background spectra. The 
resultant spectrum was grouped such that the detection 
signihcance in each data bin is at least 3cr. The confi¬ 
dence level for all the error bars reported in this paper 
are 90%. 

Another cloud feature, GO.66—0.13 

(R.A.=17^47™41.5% Decl.=-28'*26™23.0^ J2000) 


(Ponti et al. 2014), is only captured in the second 


observation (obslD 4001201901) and avoids the stray 
light from GX 3-1-1 on FPMB. We thus used both FPMA 
and FPMB from the second observation for its spectral 
analysis. The spectra were combined and grouped with 
the same method as for the central 90" region. 

2.2. XMM-Newton Data 


Instrument 

Observation 

Start Time 

Exposure 


ID 

(UTC) 

(ks) 

NuSTAR 

40012018002 

2013-10-22T16:56:07 

142.2 

NuSTAR 

40012019001 

2013-10-25T22:31:07 

151.6 

XMM-Newton 

0112971501 

2001-04-01T00:25:ll 

9.2 

XMM-Newton 

0030540101 

2002-09-09Tll:ll:26 

19.3 

XMM-Newton 

0203930101 

2004-09-04T02:53:45 

48.5 

XMM-Newton 

0694640601 

2012-09-06T10:56:15 

66.6 

XMM-Newton 

0694641301 

2012-09-16T18:34:18 

72.9 

XMM-Newton 

0694641401 

2012-09-30T19:39:50 

58.3* 


Note. — *In the observation 0694641401, the exposure time is 
58.3 ks for Sgr B2 and 44.8 ks for GO.66-0.13, as for the latter there 
are no MOSl data available. 

We collected and analyzed all the XMM-Newton 
data available in the archive covering the Sgr B2 and 
GO.66—0.13 regions. This includes observations per¬ 
formed in 2000, 2001, 2002, 2004 and 2012. The 2000 
data were excluded due to low effective exposure and 
poor statistics, while for the other years we used all the 
available observations. The list of XMM-Newton ob¬ 
servations is presented in Table 1, along with the to¬ 
tal EPIG pn-equivalent exposure times (i.e., computed 
assuming a 0.4 ratio between the effective areas of the 
MOS and pn detectors). For each selected observation 
we extracted the spectra from all available EPIC instru¬ 
ments using the XM M-Newton Extended S ource Anal¬ 
ysis Software (ESAS; |Snowden et al.||2008|) distributed 
with version 12.0.1 or the XMM-Newton Science Anal¬ 
ysis Software. For each exposure, calibrated event files 
were produced with the tasks epchain and emchain and 
filtered with pn-filter and mos-filter in order to exclude 
the time intervals affected by soft proton contamination. 
The spectra were then extracted with the ESAS mos- 
spectra and pn-spectra scripts and rebinned to have at 
least 30 counts in each bin to apply chi-square statistics. 
The quiescent component of the EPIC internal particle 
background (QPB) was estimated using archival obser¬ 
vations provided within the ESAS database and taken 
with the filter wheel closed. 

3. SPATIAL DISTRIBUTION OF THE HARD X-RAY 
EMISSION AND THE FE Ka LINE EMISSION 

Figure 1 shows the resultant 17' x 11' NuSTAR sky 
mosaics of the Sgr B2 region in the 3-40 keV, 3-10 keV, 
6.2-6.6 keV and 10-40 keV bands. The 3-40 keV image 
shows that two features are significantly detected: the 
central 90" of Sgr B2 and a newly discovered cloud fea¬ 
ture GO.66—0.13, whose 6.4 keV Fe Ka emission turned 
bright in 2012 as revealed by the XMM-Newton data 
(Terrier et al. in prep.). The green circle of 90" radius 
outlines the central region of Sgr B2, correspondi ng to 
3.4T0.3 pc with the cloud distance of 7.9±0.8 kpc ( Reid| 


et al. 2009). The green ellipse outlines the cloud feature 
GO.60—0.13, with a semi-major axis of 130" (4.9±0.5 pc) 
and a semi-minor axis of 76" (2.9±0.3 pc). GO.66—0.13 is 
^ 14 pc away from the center of Sgr B2 in the projected 
plane. Sgr B2 and GO.66—0.13 are both about 100 pc 
away from Sgr A* in the projected plane. The lower en¬ 
ergy 3-10 keV image is also shown to compare with pre¬ 
vious observations of Sgr B2 by Chandra, XMM-Newton, 
and other imaging observatories. The 6.2-6.6 keV band 
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image with continuum emission subtracted shows the Fe 
Ka line emission morphology of the Sgr B2 region. The 
10-40 keV band image provides the line-free continuum 
emission morphology. All the images are scaled indi¬ 
vidually to illustrate the morphology of major features. 
The images are overlaid with 6.4 keV line emission con¬ 
tours. The contours are made from the 2012 image of 
the XMM-Newton CMZ scan in the 6.28-6.53 keV band, 
from which the continuum emission (estimated assum¬ 
ing a power law model between two adjacent bands) has 
been subtracted (|Ponti et al.||2014| 2015 in press). 

The 10-40 keV image demonstrates that the high en¬ 
ergy X-ray (> 10 keV) morphology of Sgr B2 is re¬ 
solved at sub-arcminute scales for the first time. It 
clearly reveals substructures of Sgr B2 and proves that 
the X-ray emission in this energy band is extended (Fig¬ 
ure 1). Both the central 90" of Sgr B2 and GO.66—0.13 
are significantly detected above 10 keV. GO.66—0.13 
shows two bright cores separated by about 100", well 
correlated with its 6.4 keV line contour. The size of 
each peak is ~ 20". The 10-40 keV surface bright¬ 
ness of the whole GO.66—0.13 region is (3.5 ± 0.8) x 
10“® ph cm“^ s“^ arcmin“^. Within the central 90" 
region of Sgr B2, the 10-40 keV emission peaks at the 
center, coinciding with the compact star-forming core 
Sgr B2(M), with a surface brightness of (1.6 ± 0.1) x 
10“^ ph cm“^ s“^ arcmin“^ for the central 25" ra¬ 
dius region. It is detected at the ~ 17cr level and 
likely to be the hard X-ray counterpart of Sgr B2(M). 
The right panel of Figure 2 shows the zoomed-in T x 
5' image of the central region in 10-40 keV. Besides 
Sgr B2(M), X-ray emission around the additional com¬ 
pact core Sgr B2(N), about 50" north of Sgr B2(M), is 
also detected (~ 5cr). The zoomed-in images in Fig¬ 
ure 2 are overlaid with the source regions of Sgr B2(M) 
(RA= 17'^47™20.30*, Dec= -28°23'04.01", J2000) and 
Sgr B2(N) (RA= 17'"47™20.30^ Dec= -28°23'04.01", 
J2000 ) defined from subm m band observations by Her- 
schel (Etxaluze et al.|2013). X-ray emission from a third 
compact core Sgr B2(S), south of Sgr B2(M), is not de¬ 
tected in 10-40 keV. The surrounding regions show lower 
surface brightness, with the western half of the annulus 
from 25" to 90" brighter than the eastern half. The 10- 
40 keV morphology strongly resembles the optical depth 
map at 250 /rm, which indicates the local column density, 


deriv ed from submm continuum emission (Etxaluze et al. 
2013|). In the 250 optical depth map, Sgr B2(iVl) 


and Sgr B2(N) have highest optical depth at 250 ^m 
(T 250 Mm > 1)) while the surrounding areas have gradu¬ 
ally lower T 250 tiTn with the western region higher than 
the eastern. It suggests that the hard X-ray continuum 
emission traces the local column density of the cloud ma¬ 
terial. 

In 6.2-6.6 keV, only the main compact core Sgr B2(M) 
is significantly detected. Sgr B2(N) is not detected at 
6.4 keV, nor in the 3-10 keV range (left panel of Fig¬ 
ure 2). This could be due to higher local absorption. 
GO.66—0.13 is also not detected by NuSTAR at 6.4 keV. 
This is a dramatic change from its 2012 Fe Ka line 
morphology represented by the cyan contours, where 
GO.66—0.13 was brighter than the Sgr B2 core region. 
This indicates that the GO.66—0.13 Fe Ka emission has 
a short life time (e-folding decay time) of about 1 year. 


Finally, we note that the small white circle with 
16" radius illustrates a bright point source GXOUGG 
J174652.9—282607, which is registered i n the Chan¬ 
dra point source catalogue (|Muno et al.||2009|). The 
source is detected up to ~ 40 keV by J)JuSTAk\k.A.= 
17^47’"52.968®, Decl.= -28^26™07.37®, J2000), and will 
be reported in another work (Hong et al. in prep.). 

4. TIME VARIABILITY OF CENTRAL Sgr B2 AND 
GO.66-0.13 

With the XMM-Newton data in 2001, 2002, 2004 and 

2012 and the most recent 2013 NuSTAR observation of 
Sgr B2, we obtained the observed 6.4 keV line flux for 
the central 90" radius region of Sgr B2 using a simple 
model wabs*(apec+wabs*pow+gauss+gauss) (see Sec¬ 
tion 5.1.1) as the observed line flux is model indepen¬ 
dent. All the model parameters are fixed to be the same 
for all the years, while the normalizations for the power- 
law and Gaussian components are left free. The fit is 
satisfactory with xl = 1-0 for d.o.f. = 517. The resul¬ 
tant absorbed flux change of the central 90" of Sgr B2 
over time is shown in Figure 3 with square data points. 
The 6.4 keV line flux shows a clear decreasing trend from 
2001 to 2013, which can be fit with a linear decay model 
with a slope of s = —0.25 x 10“® ph cm“^ s“^ yr~^- The 
varying flux is preferred at a ~ 7 (t level over a flat light 
curve. In 2013 the observed 6.4 keV line flux is down 
to (0.83 ± 0.21) X 10“^ ph cm“^ s“^, which corresponds 
to about 20% of the 2001 observed flux. The 6.4 keV 
line flux decrease rate is compatible with previous works 
showing that the 6.4 keV line emission of a more ex- 
tended Sgr B2 re gion in 2005 is about 60% that of 2000 
(Inui et al. 2009). The life time of the X-ray photons 
is t ~ 11 years, which agrees well with the ~ 11-year 
light crossi ng time for the ce ntral 90", with a distance 
of 7.9 kpc (Reid et al. 2009). Although the 2001-2013 
lightcurve shows a decreasing trend, we note that the 

2013 Fe Ka emission is at the same level as that of 2012. 
With limited statistics, it is not conclusive whether the 
Sgr B2 Fe Ka emission was still decreasing in 2013 or it 
had reached a constant background emission level. Fu¬ 
ture observation of the Sgr B2 region will constrain its 
timing variability. 

Seven-year monitoring of Sgr B2 by INTEGRAL/IBIS 
reveals that the hard X-ray continuum emission flux de¬ 
creases with a similar trend to the 6.4 keV emission, 


both decay ing up to 40% from 2003 to 2009 (Terrier 
et al. 2010). This decay profile predicts that the hard 
X-ray emission in 2013 reaches ^ 30% that of 2001. 
The 10-40 keV flux of the central 90" of Sgr B2 mea¬ 
sured by NuSTAR in 2013 is Fio_4o keV = (1-9 i 0-2) x 
10“^^ erg cm“^ s“^. Therefore, the extrapolated 10- 
40 keV flux of the central Sgr B2 region in 2001 is 
Fio-40 keV 6.3 X 10“^^ erg cm“^ s“^ based on the 
hard X-ray decay profile. To verify the estimated flux 
in 2001, we extrapolated the 2001 XMM-Newton spec¬ 
trum into higher energies and derived the 10-40 keV 
flux of the central 90" in 2001 to be Fio -40 keV = 
(6 ± 2) X 10“^^ erg cm“^ s“^, which matches with the 
2001 flux value extrapolated based on NuSTAR measure¬ 
ments. This will be used to derive the luminosity of the 
primary source in Section 6.1. 

The newly discovered cloud feature GO.66—0.13, 6' 
(14 pc) from the central Sgr B2 region in the projected 
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Figure 1. Upper left panel: 17^x11' NuSTAR 3-40 keV mosaic image of the Sgr B2 region overlaid with 6.4 keV line emission contours 
(cyan dashed) made from 2012 XMM-Newton observations of the same region. The green circle is the central 90'' of Sgr B2 and the green 
ellipse shows the new cloud feature GO.66—0.13. Both the Sgr B2 core and GO.66—0.13 are detected in 3-40 keV. Upper right panel: 
NuSTAR 3-10 keV mosaic image of the Sgr B2 region.The small white circle shows a bright point source GXOUGC J174652.9—282607. 
Lower left panel: NuSTAR 6.2—6.6 keV mosaic image of the Sgr B2 core region from which the continuum emission is subtracted. The 
inner 90" of Sgr B2 correlates well with its 2012 emission. While GO.66—0.13 was brighter than the inner 90" of Sgr B2 in 2012, it is not 
detected in the 6.2-6.6 keV energy band by NuSTAR in 2013. Lower right panel: NuSTAR 10-40 keV mosaic image. The central 90" of 
Sgr B2 and GO.66—0.13 are both detected in 10-40 keV. All the images are in cts s“^ and overlaid with the Galactic coordinates with a 
grid of 0.1°. Sgr A* is outside of the field of view to the bottom right. The dark regions with irregular shapes are those of no exposure, 
which are cut out due to contamination of stray-light from bright point sources outside the field-of-view. 
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Figure 2. Zoomed-in 6'x 4' NuSTAR mosaic image of the central Sgr B2 region in 6.2-6.6 keV with continuum emission subtracted (left 
panel) and 10-40 keV (right panel). Both images are in cts s“^ and overlaid with 6.4 keV line emission contours (cyan dashed). The two 
main compact cores Sgr B2(N) and Sgr B2(M) are outlined with 15" radius region and 20" radius region, respectively. The region for the 
central 90" of Sgr B2 is overlaid for comparison. Sgr B2(M) is detected in both energy bands, while Sgr B2(N) is only detected above 
10 keV. 
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Figure 3. Time variability of the Fe Ka line for the central 90" of 
Sgr B2 and the additional cloud feature GO.66—0.13. The Sgr B2 
absorbed 6.4 keV line flux was measured in 2001, 2002, 2004 and 
2012 by XMM-Newton (red squares), and in 2013 by NuSTAR 
(orange square). The Fe Ka line flux shows a clear decay of up to 
80% of the measured Sgr B2 flux from 2001 to 2013. A chi-square 
hypothesis test favors a linearly decreasing flux at a 6.7 (t level, the 
red line shows the best-fit linearly decay model with a slope of 
s = —0.25 X 10“® ph cm“^ s“^ yr“^. The absorbed 6.4 keV line 
flux of GO.66—0.13 was measured in 2001, 2004 and 2012 by XMM- 
Newton (blue stars) and in 2013 by NuSTAR (green upper limit). 
The GO.66—0.13 Fe Ka emission experienced an increase before 
2012 and a fast decay from 2012 to 2013. The blue line shows the 
best-fit linearly increasing model with a slope of 0.12. Yet in 2013 
the 6.4 keV line was not significantly detected by NuSTAR. 

plane, was not significantly brighter than surrounding 
regions at 6.4 keV until 2012. We derived its 6.4 keV 
line flux from XMM-Newton observations in 2001, 2004 
and 2012 (stars in Figure 3). The 2002 XMM-Newton 
observation has very poor statistics for this region and 
thus was not used. The 2001 flux is poorly constrained 
to ^ 6,4 keV = 0-05 X 10“® ph cm“^ s“^ with an upper 
limit of 1.60 x 10“® ph cm“^ s“^. The 2001 and 2004 Fe 
Ka emission of GO.66—0.13 remained at the same level as 
the surrounding regions. The 2012 XMM-Newton obser¬ 
vation of Sgr B2 clearly revealed a significantly increased 
6.4 keV flux, twice that of 2004 and higher than the cen¬ 
tral Sgr B2 area. The 6.4 keV line emission contours 
shown in Figure 1 also shows that the brightest subre¬ 
gions within GO.66—0.13 were brighter than the Sgr B2 
core at 6.4 keV in 2012. The GO.66—0.13 Fe Ka emission 
experienced an increase prior to 2012 and a fast decay 
from 2012 to 2013. However, the 6.4 keV line flux sharply 
decreases and results in a non-detection by NuSTAR in 
2013, giving an upper limit (90% confidence level) to the 
6.4 keV line flux of 5 x 10“® ph cm“^ s“^. The 6.4 keV 
line flux in 2013 is less than 50% of the value measured 
by XMM-Newton in 2012. It requires a short lifetime of 
~ 1 year, which roughly matches the light crossing time 
of the two bright cores within GO.66—0.13 (2-3 years). 

5. BROADBAND 2013 SPECTRA OF CENTRAL SGR B2 
AND GO.66-0.13 

5.1. Central 9(f' Radius Region of Sgr B2 

We extracted a NuSTAR spectrum from the cen¬ 
tral 90" radius region of Sgr B2 (R.A.=17^47”"20.4^ 
Decl.=-28'‘23'"07.0", J2000). The background sub¬ 
traction method leaves mainly the molecular cloud 
emission and previously detected X-ray point sources. 


Thus, we first checked the flux contribution from 
known point sources in the source and background 
regions. The background regions we use do not 
contain any point sources above the Chandra de¬ 
tection threshold. In the source region, there are 
four point sources detected by Chandra within 90" 
of Sgr B2: CXOUGC J174723.0-282231, CXOUGC 
J174720.1-282305, CXOUGC J174718.2-282348 and 
CXOUGC J174713.7-282337. We checked each source 
using the stacked spectrum made from all archived Chan¬ 
dra data between 1999-09-21 and 2012-10-31 available 
for Sgr B2. Their summed observed flux corresponds to 
^ 7% of the total observed flux of the central 90" re¬ 
gion below 10 keV. By extrapolating their spectra into 
a higher energy band, we estimate their flux contribu¬ 
tion above 10 keV is only ^ 4%. Therefore the flux 
contribution from point sources is not significant for our 
study. None of the four sources exhibits a 6.4 keV line 
feature, while the brightest one among the four point 
sources exhibits a significant 6.7 keV line feature in its 
spectrum. It can be best fit with a collisionally ionized 
plasma model apec with a temperature of ~ 3 keV. We 
thus used this model to represent unresolved point X- 
ray sources and any possible residual GC diffuse X-ray 
background within the source region. 


5.1.1. Ad hoc XRN Model 


To examine the XRN scenario, we first applied a pop¬ 
ular XRN model, with the wabs absorption model for 
a direct comparison with previous work. The model is 
composed of a power-law and two Gaussians represent¬ 
ing the Fe Ka emission line at 6.40 keV and the K/3 
line at 7.06 keV, modified by intrinsic absorption col¬ 
umn density. All the model components are subject to 
the foreground absorption, resulting in the XSPEC model 
wabs*(apec+wabs*powerlaw+gauss+gauss) . We fixed 
the line energies of the two Gaussians to be 6.40 keV 
(best-fit centroid energy of the Fe Ka line of ^ 6.44 keV) 
and 7.06 keV, with the normali zation ratio of KfdjKa set 
to 15% (Murakami et al.||2001). 

The 3-79 keV spectrum is well-fit with this simple 
model (x^=0.97 for d.o.f. = 91). The best-fit model 
parameters are listed in Table 2. The temperature of the 
thermal component is 2 ± 1 keV with an unabsorbed 2 — 
10 keV flux of F’ 2 _io keV = (6±4) x 10“^^ erg cm“^ s“^, 
consistent with the spectra of the point sources within 
90" of Sgr B2. The 6.4 keV Fe Ka line has an ob¬ 
served flux of F 6.4 keV = (8.3 ±2.1) X 10“® ph cm“^ s“^. 
We calculated its equivalent width (EW) based on the 
power-law as the only continuum component, resulting 
in EW= 1.2lg'3 keV. The best-fit power-law photon 
index is F = 1.9 ± 0.5 , consistent w i th pr evious mea¬ 
surements of F « 2 (Terrier et al. 2010). The ob¬ 
served 10-40 keV flux is F’io- 4 o keV = (1-9 ± 0.2) x 
10“^^ erg cm“^ s“^. The intrinsic absorption column 
density is found to be A^h(*) = (5.0 ± 1.3) x 10^^ cm“^, 
on the lower end of, but still consistent with, the previ¬ 
ous result of (6.8 ± 0.5) x 10^^ cm“^ derived using the 
combined 2003 XMM-Newton and 2003 — 2004 INTE- 


GRAL/IKIS data (Terrier et al. 2010). All components 
are subject to a foreground interstellar column density of 
Nii{f) = (1.1±0.4) x 10^^ cm“^. The foreground column 
density value is consistent with our analysis of accumu- 
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lated 2001-2012 XMM-Newton data of the inner 90" of 
Sgr B2, which gives Nnii) = (1.0 ± 0.2) x 10^^ cm“^. 

Although this ad hoc model can fit well to the data, it 
is not self-consistent, with the continuum emission and 
fluorescence lines decoupled. A power-law can only mea¬ 
sure the spectral slope of the observed scattered contin¬ 
uum, but not the illuminating source spectrum. The 
model is valid for measuring the illuminating source 
spectrum only when a molecular cloud is optically thin 
{Nh ^ 10^“^ cm“^) and the Compton scattering is negli¬ 
gible. Interpretation of the resultant best-fit value for the 
intrinsic column density calls for caution. As this 

ad hoc model does not take cloud geometry into account, 
the Nnii) measured by this model represents a charac¬ 
teristic column density of the cloud, while in reality il¬ 
luminating X-ray photons are absorbed and scattered in 
various locations of the cloud. 


5.1.2. Self-consistent XRN Model MYTorus 

In the XRN scenario, to consistently measure the illu¬ 
minating X-ray spectrum and to properly determine the 
intrinsic column density, a self-consistent XRN model 
based on Monte-Car lo simulations is required. 
rakami et al. _( |200f ), Revnivtsev et al. (2004), 


Mu- 


Ter- 


rier et al. 'i 2010p and Udaka et al. 12011) have applied 


Monte-Carlo based XRN models to Sgr B2 data in order 
to study its morphology and spectrum. The MYTorus 
model is the only XRN model available in XSPEC to self- 
consistently measure the illuminati ng source spectrum 
and the intrinsic column density (Murphy & Yaqoob 
2009). MYTorus was originally developed for Compton- 
thick Active Galactic Nuclei (AGN) assuming a toroidal 
reflector with neutral materials and uniform density. The 
default Fe abundance in the MYTorus model is one so¬ 
lar and does not allow variation. There are three com¬ 
ponents offered in the model: the transmitted contin¬ 
uum (MYTZ), the scattered continuum (MYTS) and the 
iron fluorescence lines (MYTL). In the case of GCMCs, 
the observed spectrum only contains the last two com¬ 
ponents, because we are seeing only the reflected X- 
ray photons off the cloud. We thus use the com¬ 
bination of MYTS and MYTL. Both components de¬ 
pend on the intrinsic equatorial hydrogen column den¬ 
sity Nnii), the illuminating source spectrum photon in¬ 
dex T, the inclination angle 9obs between the line-of-sight 
and the torus symmetry axis, and the model normal¬ 
ization Nmt- To self-consistently measure the illumi¬ 
nating source spectrum, all the parameters are linked 
between MYTS and MYTL as a coupled mode where 
the same incident X-ray spectrum is input into both 
components. We select a termination energy of the 
incident power-law as 500 keV. As the best-fit energy 
centroid for the Fe Ka line is ~ 6.44 keV, we select 
an energy offset of -1-40 eV for the MYTL component 
to allow freedom for the centroid energy. The resul¬ 
tant model is wabs* (apec+MYTS+MYTL) , where the two 
MYTorus components are implemented as table models 
in Xspec: atable{mytorus_scatteredH500_v00.fits}+ 
atable{mytl_V000010pEp040H500_v00.fits}. 

The coupled mode of the MYTorus model can be ap¬ 
plied to the GCMC X-ray reflection spectra with some 
restrictions, by treating a quasi-spherical cloud as part of 
a virtual torus and rescaling incident X-ray flux properly 
(see Figure A6 for geometry). Due to the toroidal geom¬ 


etry and the uniform density profile MYTorus assumes, 
for GGMGs we restrict the applicability of the model to 
the inclination angle range of dobs ^ 60° and the equato¬ 
rial column density range of Nh ^ 10^^ cm“^, in order 
to be insensitive to the torus geometry (see Appendix 
for more details). When Nnii) reaches ^ 10^"^ cm“^, we 
derive the systematic error based on the MYTorus model 
to be 25% for Nh, 3% for F and 10% for the model nor¬ 
malization (Mori et al. submitted). We note that the 
Nnii) measured by the MYTorus model is an averaged 
value over the torus. 

This self-consistent model can fit to the Sgr B2 spec¬ 
trum well, yielding xi = 1-09 for d.o.f. = 92 (left 
panel of Figure 4). Both the best-fit temperature and 
the flux of the apec component are consistent with 
the fitting results of the ad-hoc power-law model (Ta¬ 
ble 2). The foreground column density is Nnif) = 
(1.2 ± 0.1) X 10^^ cm“^, consistent with that derived 
from the ad hoc model and better constrained. How¬ 
ever, the intrinsic equatorial column density is Nnii) = 
(1.01 ± 0.16stat ± 0.25sys) X 10^^ cm“^, which is twice 
the Nnii) value measured with the ad hoc model. The 
Nnii) value derived with the two models cannot be di¬ 
rectly compared due to the lack of geometry definition 
for Nnii) in the ad hoc model. Ab^(^) in the MYTorus 
model corresponds to the minor diameter of the torus 
(or diameter of the quasi-spherical cloud), while Nnii) 
in the ad hoc model roughly corresponds to the cloud ra¬ 
dius as it “measures” an averaged effect, which might 
explain the difference by a factor of 2 (see the Ap¬ 
pendix). For the central 90" of Sgr B2, the cloud ra¬ 
dial optical depth due to Thomson scattering derived by 
the MYTorus model is tt = 0.67 ± 0.27, consistent with 
previous measurements of ry = 0.4 ± 0.1 derived from 
a self-developed Monte-Garlo simulation with combined 


ASCA/GIS, GRANAT/AYd-V and INTEGRAL/imS, 


data when Sgr B2 was much brighter (Revnivtsev et al., 
2004|). As a result, the central 90" of Sgr B2 is marginally' 


optically thin. The illuminating source spectrum photon 
index is constrained to be F = 2.2 ± 0.4 with the sys¬ 
tematic error negligible, compatible with the observed 
spectrum slope of F = 1.9 ± 0.5 measured by the power- 
law model which does not consider Compton scatter¬ 
ing. With the hard upper limit on the inclination an¬ 
gle Oobs = 60° chosen for this study, the MYTorus model 
is not sensitive to the parameter Oobs- Therefore, this 
parameter could not be constrained by our work. 

Though with limitations, MYTorus is currently the best 
self-consistent XRN model available. The systematic er¬ 
rors could be underestimated since they are estimated 
solely based on the MYTorus model (see the Appendix). 
A modified XRN model based on MYTorus with a spheri¬ 
cal reflector and variable Fe abundance is under develop¬ 
ment for the specific case of GCMCs, which will remove 
the uncertainties induced by the toroidal reflector and 
the fixed Fe abundance (T. Yaqoob, private communi¬ 
cation). The more sophisticated XRN model will also 
implement the specific density profile of Sgr B2 to re¬ 
move uncertainties induced by non-uniform density and 
better constrain Nnii) (Walls et al. in prep.). 


5.1.3. Self-consistent Cosmic Ray Electron Model LECRe 
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To determine whether CR electrons can be a ma¬ 
jor contributor to the remaining Sgr B 2 X-ray emission 
in 20 13, we tested the LECRe model (Tatischeff et al. 
20121. The LECRe model has four free parameters: the 
power-law index of the accelerated particle source spec¬ 
trum, s, the path length of cosmic rays in the non- 
thermal X-ray production region, A, the minimum elec¬ 
tron energy, Emin, and the model normalization, Nlecr, 
which provides the power injected in the interaction re¬ 
gion by primary cosmic-ray electrons of kinetic ener¬ 
gies between Emin and Emax = 1 GeV: dW/dt = 
^ttD^Nlecr erg s“^, where D is the source distance. 
Since the fit cannot constrain the path length of LECRe, 
we fix it to A = 5 X H — atoms cm“^ as was used in 


Tatischeff et al. 


2012 

E'lgure 4, it resu 

ts in 


Xy = 1.18 for d.o.f. = 92 (Table 2). The electron spectral 
index is found to be s = 2.8±0.5 with a normalization 
of Neecrc = (2.8 ± 0.2) X 10“^ erg cm“^ s“^. Previ¬ 
ous analysis with XMM-Newton and INTEGRAL/ISIS 
data found that s 1.5 , but with an unconstrained error 
bar (Terrier et al.||2010). The best-fit minimum electron 
energy is as low as 3 keV, with an upper bound of 50 
keV. Electrons with such low energy cannot leave their 
acceleration site, nor penetrate the cloud. The best-fit 
molecular cloud metallicity is Z/Zq = d.O+g^e’ higher 


than ZjZQ = 3.1 ± 0.2 derived by Terrier et al. (2010). 
Such a metallicity value is too high compared to current 
measurements of the GC interstellar medium metallic¬ 
ity which ranges from slightly higher than solar to twice 
solar. Furthermore, compared to the XRN models, the 
LECRe model results in a poorer fit above 10 keV as it 
cannot fit to the spectrum curvature equally well. The 
high energy part of the spectrum thus provide a more ex¬ 
cluding constraint on the LECRe model, which does not 
depend on the high metallicity required by the signifi¬ 
cant Fe Ka line. We therefore confirm that the unphys¬ 
ical best-fit model parameters makes the LECRe scenario 
unlikely to be a major process to account for the Sgr 
B2 X-ray emission in 2013. If there is an underlying con¬ 
tribution from CR electrons, it has to be significantly 
lower than the current level. 

5.1.4. Self-consistent Cosmic Ray Proton Model LECRp 

While CR electrons can be safely excluded as a ma¬ 
jor contributor to the 2013 Sgr B2 X-ray emission, low 
energy CR proton/ion bombardment could be a major 
process if the Sgr B2 X-ray emission obtained in 2013 
has already reached the constant background level. The 
LECRp model has the same model parameters as the 
LECRe model, with four free parameters s. A, Emin and 
Nrecr- As the path length cannot be constrained by the 
fit, we fix it to A = 5 X 10^"^ H — atoms cm“^, a typical 
value for nonrelativistic particles propag ating in massive 
molec ular clouds of the GC environment (Tatischeff et al. 
2012). The right panel in Figure 4 shows that it results 
in an overall good fit with = 1.17 for d.o.f. = 92 
(Table 2). The best-fit Fe abundance is Z/Zq = 2.51j;q, 
which is consistent with the GC metallicity. The CR 
proton spectral index is s = 1.9^0®. This agrees with 
the CR proton spectral index derived with the LECRp 
model for the arch es region using the XM M-Newton spec¬ 
trum (s = 1.9^Q g, Tatischeff et al. [2012 ) or the NuSTAR 


spectrum (s = 1.7 ± 0.6, Krivonos et al. |2014 ) CR iron 
bombardment was proposed by these autnors to explain 
the X-ray emission from the Arches before a sig nificant 
variability of the X-ray e mission was detected ( Clavelf 
et al. 120141). As discussed in|Tatischeff et al.|(|2012 ), with 
such a relatively hard CR spectrum, the X-ray emission 
produced by CR protons depends weakly on the mini¬ 
mum ion energy Emin- Therefore, we fix it to Emin = 
10 MeV nucleon”^. The total power required by CR pro¬ 
tons in the cloud can be derived from the best-fit model 
normalization Neecr = (1.4±0.4) x 10“^ erg cm“^ s“^. 
The power injected by primary protons of energies be¬ 
tween Emin = 10 MeV and Emax = 1 GeV is dW/dt = 
(1.0 ± 0.3) X 10®® erg s“® for D = 7.9 keV. However, 
the power injected into the cloud by the CR protons de¬ 
pends on the minimum energy Emin- For Emin = 1 MeV 
and Emin = 100 MeV, the best-fit model normaliza- 

4-11 10-7 QT^-2 g-l 

-7 „^-2 „-i 


tions are Neecr = 2 . 0 _g ;5 x 10 


Neecr = 0.7_qX x 10“' erg cm“^ s“®, respectively. 
The corresponding power injected by CR protons is in 
the range of dW/dt = (0.4 — 2.3) x 10®® erg s“®. There¬ 
fore, the LECRp scenario could explain the current level 
Sgr B2 X-ray emission and can only be excluded by fur¬ 
ther variability. 

5.2. Main Compact Cores Sgr B2(M) and Sgr B2(N) 

Herschel observations have shown that the local 250^m 
optical depth at Sgr B2(M) and Sgr B2(N) (T 250 ;im ~ 1) 
is higher than that of the s urrounding region w ithin 90" 
roughly by a factor of 2-5 (Etxaluze et al. ]M^ . Specif- 
icall y, the Sgr B2(M) 250 /im optical depth is T 2 r,viu.m ^ 
0.9. Etxaluze et al. (2013) converted it to a local column 
density of Nh ~ 1 x 10" cm“® based on the the conver¬ 
sion relationship T 250 tJ.m/NH = 8.8 x 10“®® cm® for opti¬ 
cally thin clouds. This conversion factor is increased to 
5 X 10“®® cm® for the Sagittarius region (Bernard et al. 


2010), thus reducing the Sgr B2(M) column density to 
Nh ~ 1.8 X 10®'^ cm“®. We extracted a spectrum from 
a 20" radius region (Figure 2) centered on Sgr B2(M) 
to measure the local optical depth independently using 
the X-ray data. The Sgr B2(M) spectrum is fit with the 
MYTorus model using the same model settings for the 
central 90" region. The model yields a good fit with 
X® = 1.06 for d.o.f. = 39. While all the other key pa¬ 
rameters are consistent with those derived from the cen¬ 
tral 90" region, the intrinsic equatorial column density 
is Nnii) = (2.0('(^(g()“)('(Q g(“() x 10®“^ cm“®, twice that of 
the whole 90" region but compatible within error bars. 
The LECRP model can also fit well to the Sgr B2(M) spec¬ 
trum (x® = 1.02 for d.o.f. = 39), giving consistent Nnii) 
value. It agrees well with the Herschel measurements of 
local column densities. The corresponding radial opti¬ 
cal depth in the direction of Sgr B2(M) is tt = 1.3tj;g. 
The NuSTAR observations confirm the Herschel mea¬ 
surements that the Sgr B2(M) region is not optically 
thin, though with large error bars. The 10-40 keV flux 
of Sgr B2(M) and Sgr B2(N) falling into the central 90" 
region is F’ig_ 4 g = (6.1 ± 0.2) x 10“®® erg cm“® s“®, 
contributing to about one third of the total flux of the 
central 90". The surrounding regions within the 90" ra¬ 
dius region contribute to the remaining two thirds of the 
total flux. The second core, Sgr B2(N), does not have 
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sufficient statistics to perform a similar analysis. 

5.3. The new feature GO. 66—0.13 

Since neither focal plane modules are severely con¬ 
taminated by stray light, we were able to use the com¬ 
bined FPMA and FPMB NuSTAR data for GO.66-0.13. 


et al. 

2001 

Nowak et al. 

2012 

Barriers et al. 

2014 

), al 

though the mechanisms o 

both the past giant 

X.-ray out 


bursts and the current relatively mild X-ray flares remain 
unclear. The outburst spectrum is also compatible with, 
though at the extreme end of, the average AGN photon 


Furthermore, the GO.66—0.13 region does not contain 

to 300 keV 

(|Molina et al. 

2013 

), which is also expected 

sources registered in the Ghandra point source catalogue 

for the Sgr 

4* spectrum iJ 

it became active. 


(Muno et al. 2009). The background-subtracted NuS- 


The Sgr B2 cloud as a whole has been discussed as an 


lAR spectrum tor GO.66—0.13 does not show any signif- 

optically thin cloud (e.g. |Koyama et al. 1996 Sunyaev 

leant line features. To compare with the significant Fe 

& ChurazovI 

1998[ Revmvtsev et al.||20U4[ I'Ferrier et al. 


Kof lines detected by XMM-Newton in 2012, we fit the 
NuSTAR 2013 spectrum together with the XMM-Newton 
2012 spectrum. Since the sharply decreasing Fe Ka emis¬ 
sion cannot be explained by the LEGRp scenario, we ex¬ 
amine the GO.66—0.13 spectrum only with the MYTorus 
model. Because of different background subtraction 
methods, the thermal component in the XMM-Newton 
data represents the GG diffuse emission while the ther¬ 
mal component in the NuSTAR data represents unde¬ 
tected point source and residual diffuse emission. There¬ 
fore we use the model wabs* (apec+MYTS+MYTL) for the 
joint fitting, with all the model parameters linked except 
for kT and the model normalizations of apec, MYTS and 
MYTL. The model gives a good fit (x^=1.18 for d.o.f. = 
55), with a best-fit photon index of F = 1.4± 0.5, a fore¬ 
ground column density of Nnif) = ( 8 . 2 j'j 4 ; 5 ) x 10 ^^ cm“^ 
and an intrinsic equatorial column density of = 

3.0l4;g X 10^^ cm“^ (Figure 5). The intrinsic density 
Nffii) is lower than that of the central 90" radius re¬ 
gion, indicating that GO.66—0.13 is optically thin. How¬ 
ever, in case that GO.66—0.13 is partially illuminated, 
the measured which is the illuminated column 

density, would be smaller than the cloud intrinsic col¬ 
umn density. The 10-40 keV flux for the NuSTAR 2013 
data is Eio_ 4 o keV = (9-0 i 1 - 2 ) x 10 “^^ erg cm“^ s“^ 
in 10-40 keV. The observed 6.4 keV flux in 2012 mea¬ 
sured by XMM-Newton is Fqa keV = (1-5 ± 0.2) x 
10“^ ph cm“^ s“^, while the 2013 NuSTAR spectrum 
does not require a 6.4 keV line, with 6.4 keV flux upper 
limit of 5x 10“® ph cm“^ s“^ (90% confidence level). The 
8-12 keV non-thermal continuum emission measured by 
NuSTAR in 2013 dropped to ^ 50% of that measured 
by XMM-Newton in 2012. Both the fluorescence emis¬ 
sion and the hard X-ray emission show fast variability 
(or short life time), which could be due to a partial illu¬ 
mination of GO.66—0.13. 

6. DISCUSSION 

6 . 1 . X-ray Reflection of Past Sgr A* Outbursts 

The broadband spectrum obtained by NuSTAR sup¬ 
ports the XRN scenario where Sgr B2 reflects incom¬ 
ing X-rays from an illuminating source, which is most 
likely a past X-ray outburst from Sgr A*. The illuminat¬ 
ing source photon index measured with the broadband 
NuSTAR spectrum in 2013 using the MyTorus model is 
constrained to F = 2.2 ±0.4, compatible with the illumi¬ 
nating source photon index of F = 2.0 ± 0.2 derived with 
combined XMM-Newton and INTEGRAL/IBIS spectra 
obtained in 2003-2004 when Sgr B2 was much more 
brighter (Terrier et al. 2010). Such a past Sgr A* out¬ 


burst spectrum is consistent with the curre nt bright 
Sgr A* X-ray flares spectra with F « 2 —3 (e.g. Baganoff 


2010[ ). in 2013, NuSTAR found the harci X-ray emis¬ 


sion is concentrated in the central 90" radius region of 
Sgr B2, but it is unknown whether the hard X-ray emis¬ 
sion in 2003-2009 was similarly concentrated within 90" 
or more extended due to the mild angular resolution of 
INTEGRAL/ISIS. While in the direction of the compact 
core Sgr B2(M) the optical depth is as high as tj’ > 1, 
the majority of the central 90" radius region is optically 
thin with tt < 1. Besides, the 10—40 keV image shows 
that the X-ray continuum emission traces the local col¬ 
umn density, also suggesting the majority of the central 
90" of Sgr B2 is optically thin. Therefore, although there 
are complicated substructures with higher local column 
densities, the central 90" can be considered as optically 
thin for a simplified calculation of Sgr A* outburst lumi¬ 
nosity. 

For an optically thin cloud, the luminosity of the pri¬ 
mary source can be derived via two independent analyt¬ 
ical methods: the Fe Ka fluorescence process using the 
6.4 keV line flux we measured, and the Gompton scat¬ 
tering process with the measured hard X-ray continuum. 
In most previous works, the primary source luminosity 
has been calculated via the fluorescence process, due to 
good measurements of the 6.4 keV line but poorly con¬ 
strained continuum emission. The 6.4 keV line emission 
flux for continuu m radiation is expressed by |Sunyaev fc| 
Ghurazov (11998) as 


^6.4 = 


n 


AttD^ 


ZpeTTl{8) ph cm 


-2 s-i 


( 1 ) 


where H is the cloud solid angle subtended to the illumi¬ 
nating source, D is the distance to the observer, Zpe is 
the iron abundance relative to solar, Tp is optical depth 
due to Thomson scattering, and I(8) is the illuminat¬ 
ing source flux at 8 keV. (j> is a factor of the order of 
unity, which depen ds on the source spectrum. [Sunyaev] 
& Ghurazov (1998) assumed the source spectrum to be 
bremsstrahlung with a temperature from 5 to 150 keV, 
corresponding to (/) ~ 1.0 — 1.3. Now, with a better 
knowledge of the primary source spectrum, we know it 
can be well described with a power-law with F « 2. F 
ranges from 1.8 to 2.6, as measured with the MYTorus 
model, corresponding to ^ ~ 1.1 — 1.3. In the calcula¬ 
tion hereafter we will adopt F « 2 and its corresponding 
6 = 1.18. 


Using the Sunyaev & Ghurazov (1998) approach, we 
express I(8) with illuminating source luminosity at 8 keV 
in the 8 keV wide energy band as ±8 = d( 8 ) x 8 x 8 x 
1.6 X 10“® erg s“^. For the central region of Sgr B2 we 
adopted r = 90" assuming it was entirely illuminated, 
so the luminosity of the source required to produce the 
observed 6.4 keV line is: 
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Figure 4. The 3-79 keV X-ray spectrum of the central 90^^ of Sgr B2 measured with NuSTAR FPMA fitted with three different models. 
The crosses show the data points with Icr error bars. The solid lines are the best-fit models. The lower part of the plots shows the deviation 
from the model in units of standard deviation. The spectrum is fitted with the self-consistent XRN model MYTorus (left panel), the LECRe 
model (middle panel) and the LECRp model (right panel). While all the models can fit the spectrum well overall, the LECRe model results 
in a poorer fitting above 10 keV. The LECRe scenario is ruled out based on unphysical best-fit model parameters. 



5 10 20 

Energy (keV) 


Figure 5. 2013 NuSTAR (black) and 2012 XMM-Newton (red for PN and green for MOS) spectra of GO.66—0.13 jointly fit to the XRN 
model wabs* (apec+MYTS+MYTL). The crosses show the data points with Icr error bars, and the solid lines show the best model. The lower 
panel shows the deviation from the model in units of standard deviation. The 6.4 keV line was significantly detected by XMM-Newton in 
2012, but not detected by NuSTAR in 2013. It suggests a fast decay of fluorescence emission within one year. 


Lr = 3x10 


39 




6.4 


10 


-4 


0.1 


tt ) \100 pc 


(Zpe) ^ ergs ^ 

( 2 ) 

where d is the distance between the molecular cloud and 
the illuminating source. For a power-law spectrum with 
r « 2, Lg corresponds to 30% of the source luminos¬ 
ity at 3 — 79 keV. While the 6.4 keV line flux has been 
decreasing over the past decade, the maximum observed 
Fe Ka flux should be used to calculate the peak pri¬ 
mary source luminosity. As the XRN model we use does 
not measure the iron abundance of Sgr B2, we u sed the 
most recen tly measured value Zp^. = 1.3 ± 0.1 (|Terrier 
et al.|2010 ). With the maximum central 90" Sgr B2 hux 
-^6.4 = (4.3 ± 0.6) X 10“® ph cm“^ s“^ obtained in 2001 
and averaged optical depth Tp = 0.67 ± 0.27 measured 
with by NuSTAR, the resultant primary source luminos¬ 
ity is Lg = (1.5 ± 0.7) X 10^®(d/100 pc)^ erg s“^, corre¬ 
sponding to a 3-79 keV luminosity of A3-79 keV = (5.0 ± 
2.3) X 10^®(d/100 pc)^ erg s“^ or 1-100 keV luminosity 
of Ai_ioo keV = (0.7 ± 0.3) X 10^®(d/100 pc)^ erg s“^. 
This is comparable to the Li_ioo keV ~ 1.1 x 

10^® (d/100 pc)® erg s“^ derived from XMM-New ton and 


INTEGRA L/IBIS data obtained in 2003-2004 (Terrier 
et al.|2010| . The uncertainty of the primary source lumi- 


nosity mainly comes from the measurement uncertainty 


of Tp and the distance between Sgr B2 and the primary 
source, d. The solid angle could also be lower in the 
case of partial illumination, which would result in an un¬ 
derestimation of the primary source luminosity. 

Since this was the first time the Sgr B2 10-40 keV 
continuum emission has been resolved, we can also de¬ 
rive the primary source continuum luminosity from the 
Comp t on sca ttered continuum. As discussed in |Capelh] 


et al. (2012), the observed scattered continuum at an 
angle 0 with respect to the incident radiation direction 


is: 


S(E) 


n 3 tp (1 -I- cos® 9) 
2 


I{E) ph cm-® 


keV-^ 


( 3 ) 

where I{E) is the photon flux from the source in units 
of ph s”^ keV”^. This relationship is based on Thomson 
scattering, which is a good approximation for incident 
photons with energies of a few keV. However, when the 
photon energy is higher, electron quantum mechanical 
effects become important. As a result, we need to con¬ 
sider the Compton scattering process. Starting with the 
Compton scattering cross section formula, we derived a 
correction factor from Thomson scattering to Compton 
scattering as fpc = (1 + snAv (^ ~ cosd))”"^ assuming 
a power-law spectrum with F « 2. The correction holds 
as long as hu> <C me® for the electrons, which is true for 
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Table 2 

Spectral analysis of central 90^^ of Sgr B2 with NuSTAR data. 


Parameter 

Unit 

Power-law ModeP 

MyTorus Model*’ 

LECRe ModeP 

LECRp Model** 

Nuif) 

10^3 cm-3 

1.1 ± 0.4 

1.2 ±0.1 

1.3 ± 0.5 

1.2 ± 0.6 

Nuii) 

10^3 cm-3 

5.0± 1.3 

10.1 ± 4.1 

5.5 ± 1.3 

'^•^-2.6 

ZIZq 



1(fixed) 


9 c + l-S 
^•^- 1.0 

kT 

keV 

2 1 

^•-*--0.9 

0 r-|-2.5 

o.O_i_4 


^•^-0.7 

Fapec(2 - 10 keV) 

erg cm~^ 

etl X 10-33 

fits X 10-33 

7/3 X 10-33 


tb.4keV 

10 “® ph cm“3 g-i 

8.3 ± 2.1 




EW6.4 keV 

keV 

-*-•^-0.3 




rrij. 


1.9± 0.5 




> 

0 

1 

0 

erg cm~^ 

(1.9 ±0.2) X 10-32 




^ MT 



2.2 ±0.4 



Nmt 

10“3 ph cm“3 g-i 





A 

H — atoms cm~^ 



5 X 1024 (fixed) 

5 X 102*3 (fixed) 

s 




2 -8+0 5 

1 q+0-8 

Fj 

■‘-'min 

keV 



3+47 

10^ (fixed) 

Nlecr 

10 “^ erg cm~^ s~^ 



2.8 ±0.2 

1.4 ± 0.4 

xl (d.o.f.) 


0.97 (91) 

1.09 (92) 

1.18 (92) 

1.17 (92) 


Note. — The goodness of fit is evaluated by Xu number of degrees of freedom given in parentheses. The errors are 90% 

confidence. 

“ Power-law Model: wabs*(apec+wabs*pow+gauss+gauss). 

Power-law model is characterized by foreground interstellar absorption column density 
temperature kT and 2—10 keV flux Fapec(2 — 10 keV) of apec, intrinsic absorption column density 
photon index Pp./. and 10 — 40 keV observed flux Fpj.(10 — 40 keV) of power-law, equivalent 
width EW and observed flux -F6.4keV of fh® 6.4 keV iron line. 

^ MYTorus model: wabs* (apec+MYTorus). 

The MYTorus model is characterized by A^h(/)? ^i^d Fapec(2 — 10 keV) of apec, the intrinsic 
absorption column density Nn (i) , photon index Pmt of f ho illuminating power-law spectrum 
and the power-law normalization A^mTj while the inclination angle is fixed to = 0°. 

^ LECRe model: wabs* (apec+wabs*LECRe). 

The LECRe model is characterized by Nuif), f^T and Fapec(2 — 10 keV) of apec, the power-law 

index s of the accelerated particle source spectrum varying from 1.5 to 5, the path length A of the 
cosmic rays in the non-thermal X-ray production region varying from 10^^ to 10^® H — atoms cm“^, 
minimum energy Emin and the model normalization Nlecr foi" LECRe. Nlecr provides the power 
injected in the interaction region by primary cosmic-ray electrons or protons of kinetic energies between 
Emin and Emax — 1 GeV. 

^ LECRp model: wabs* (apec+wabs*LECRp). 

The LECRp model is characterized by Nuif), ^nd Fapec(2 — 10 keV) of apec, s ranging 

from 1 to 5, A ranging from 10^^ to 10^^ H — atoms cm“^, Emin and the model Nlecr for LECRp. 

the NuSTAR energy band of 3 to 79 keV. 

Thus, for any X-ray photon energy, the Compton scat¬ 
tered continuum can be expressed as: 


(1 — cos 0)j 

(4) 

The observed 10-40 keV scattered continuum flux is 
-^10-40 = Iio keV dE erg cm”^ g-b Thus, in¬ 

tegrating both sides of equation (4) gives us the rela¬ 
tionship between Fio-40 and Lg. We keep the first or¬ 
der correction term of the integrated form and neglect 
higher order terms. As a result, the source luminosity 
in the 8 keV wide band can be expressed with observed 
10-40 keV Compton scattered continuum measured at 
angle 0 as: 


S{E) = 


n 3 tt (1 -I- cos^ 9) 
2 

xI{E) ph cm 


1 + 


E 


-2 g-i keV 


511 keV 

-1 


L8 = 4 X 10^® 

1 

X 


T'ln- 


10-40 


10“i^erg cm“2 g-i 


0.1 

tt 


f ^ 


1 -f cos^ 0 


\100 pc^ 

X (0.72-^0.12(1-cos6»)) erg s-i (5) 


We use the maximum value of the hard X-ray contin¬ 
uum flux in 2012 from the lightcurve obtained by XMM- 
Newton. In Section 5 we estimated the 10-40 keV flux 
of the central 90" in 2001. The resultant source luminos¬ 
ity is Ls = (3.6 ± 1.8) X 10^8 f{0) (d/100 pc)^ erg s-\ 
where f{0) = i+^(0.72 + 0.12(1 - cosd)) ~ 0.36- 
0.84 corresponding to 0 ~ 0°-180°. As the scattering 
angle range cannot be derived from the XRN model fit¬ 
ting in our study, we adopt cos 6> ^ 0 .8 based on VLBI 
observation results (Reid et al. 20091 to derive a nomi¬ 
nal value for f{9). The error of f[0) is calculated cor¬ 
responding the full allowed range of the scattering an- 





















12 


S. Zhang et al. 


gle e = 0°-180°, resulting in f{e) = 0.46t^;?^. Thus, 
Lg = ^ (d /100 pc)^ erg s“^, correspond¬ 

ing to a 3-79 keV luminosity of ^ 3-79 keV = (S-O^^ g) x 
10^®(d/100 pc)^ erg s“^ or a 1-100 keV luminosity of 
Ti -100 keV = (0-8to;4) X 103®(d/100 pc)^ erg s“k 

This result is consistent with Lg calculated via the pho¬ 
toelectric absorption process. Indeed, comparing equa¬ 
tion (2) and (5), we have 


Lg(a6s) 

Lg(scot) 


-f^lO-40 


10 - 


-^erg cm“ 


—) {Zpe 


i-1 


TVi 


f[e) 




Y 10“^ph cm~^ £ 

( 6 ) 

which is a factor of the order of unity. The pri¬ 
mary source luminosity independently calculated via the 
Compton scattering is remarkably consistent with that 
calculated via photoelectric absorption. The consistency 
strongly supports the X-ray reflection scenario. With 
Equation (4), the primary source luminosity in any en¬ 
ergy band can be directly derived given the measurement 
of the molecular cloud flux in the same energy band. This 
continuum-flux-based method has also been applied to 
other GC molecular clouds to test the XRN model and 
constrain the Sgr A* outburst history (Mori et al. 2015). 

The MYTorus model can also give the illuminating 
source spectrum, although it assumes fixed solar abun¬ 
dance for iron and that the observation angle cannot 
be constrained with the data. By rescaling the incom¬ 
ing X-ray flux by the solid angle ratio of the torus 
(0/47r = 0.5) and Sgr B2 {d/An = 2.89 x 10-^), 
we derive the primary source luminosity to be Lg = 
(1.4+9;®) ^ (d/100 pc)^ erg s“^ (see the Appendix 

for details). As the 2013 X-ray emission from Sgr B2 
is 30% that of its maximum flux (Section 5), the pri¬ 
mary source luminosity with the peak Sgr B2 flux is 
Lg = (d.etg;^) X 10®® (d/100 pc)^ erg s“®, or L 3-79 keV = 
(1.5^?;^) X 10®9(d/100 pc)2 erg s"! with the 2001 Sgr B2 
X-ray flux. With large error bars, the primary source 
luminosity derived by the MYTorus model overlaps with 
that derived with the analytical calculations, though the 
best-fit luminosity value is higher. For the other GC 
molecular clouds, the MYTorus model also gives a primary 
source luminosity higher than that derived from the an¬ 
alytical calculation roughly by a factor of 2 (Mori et al. 
2015). 

We note that the analytical calculations above are valid 
for an optically thin (tt 'C 1) cloud where the fluo¬ 
rescence photons are not scattered and the continuum 
emission photons are scattered once inside the cloud. 
Sgr B2 has very complicated density profile, including 
three main compact cores (Sgr B2(N), Sgr B2(M) and 
Sgr B2(S) from north to south) and numerous clumpy 
regions. Overall its density decreases as radius becomes 
larger. Without a model fully considering the density dis¬ 
tribution within Sgr B2, the measured tt for a specific 
region is an averaged value. For the 90" radius region 
we adopt, and larger regions previous wo rks adopt (e.g. 
4.5' radius region in Terrier et al. 2010), the averaged 
optical depth for such large regions is ry 1, thus have 
been treated as optically thin as a whole. However, for 
the compact cores at the center of Sgr B2, which have 


higher local density, the local optical depths are higher 
than 1 as indicated by Herschel and NuSTAR measure¬ 
ments, though with large error bars. For these regions, 
multi-scattering could take place during the decay stage. 
The X-ray emission life time of optically thick material 
is longer than its light crossing time {t = TTr/c when 
r 1) due to decay caused by multi-scattering. The 
reflected X-ray emission follows an exponential decay for 
optically thick materials and a parabolic decrease for op¬ 
tically thin materials. Therefore, the optically thicker 
material flux decreases much slower at late stages. If the 
compact core regions are confirmed to be optically thick. 
X-ray emission from the central compact cores Sgr B2(M) 
and Sgr B2(N) could still be bright enough to be detected 
when the surrounding optically-thin regions completely 
fades. However, if the optical depth in the direction of 
the compact cores can be more tightly constrained to 
Tt = 1, the multi-scattering effect would be negligible. 
Future observation of Sgr B2 can test this prediction and 
better constrain the local optical depth in the direction 
of the compact cores. 


6.2. Cosmic Ray Electron Bombardment 

LECRe cannot be the major contribution to the bright 
Fe Ka line emission from Sgr B2, as both the required 
iron abundance and t he electron power are too high 
to be phys ical (e.g. Revnivtsev et al. 2004 Terrier 
et al. 2010). After a decade of decreasing m the Fe 
Ka line flux, we tested whether the the constraints on 
the LEGRe parameters by the remaining emission level 
in 2013 are still excluding this scenario. Generating 
the bremsstrahlung emission with the observed slope of 
r = 1.9 in 3-79 keV requires a relatively soft electron 
spectra with s ^ 2.8 and a minimum GR electron en¬ 
ergy far below 100 keV as we m easured, which is illus¬ 
trated in Tatischeff et al.l (2012). For such a soft GR 
spectrum, the neutral Fe Ka line is predicted to be rel- 
atively weak, with EW < 0.4 x (Z/Z©) keV (Fig. 3 in 
Tatischeff et al. 2012|). Therefore, to fit the observed 
Fe Ka line FW of ~ 1.3 keV as measured by NuSTAR^ 
the required Fe abundance is found to be Zpe = 4.0+9(9, 
still significantly exceeding the current measurement of 
GG metallicity of solar to twice solar. This is due to 
the relatively inefficient production of 6.4 keV photons 
by LECRe interactions, for which the iron fluorescence 
yield R^a keV = Tx(6.4 keV)/(dWe/dt) is always lower 
than 1 X \t)~^(ZIZcA) for s > 2.3 and Ej^in < 100 keV 
(Tatischeff et al.|2012|. Besides the unphysical Fe abun- 
dance required by the Fe Ka emission, the LECRe model 
also results in an overall poorer fitting to the broadband 
3-79 keV spectrum compared to the XRN model. Es¬ 
pecially, the hard X-ray part of the spectrum has more 
residuals for the LECRe model, thus providing a more ex¬ 
clusive constraint which does not depend on local metal¬ 
licity. 

Another difficulty in explaining all the 6.4 keV emis¬ 
sion with LEC Re is the lack of an obvi ous particle ac¬ 
celeration site. Yusef-Zadeh et al. (2002) suggested that 
the interaction ot non-thermal radio hlaments and the 
cloud could locally accelerate GR electrons. Yet until 
now no such filaments have been observed around Sgr B2. 
Even if there exists a local GR acceleration site, the non- 
thermal electrons with energies below 100 keV are not 
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able to escape the acceleration region and p enetrate the 
whole dense cloud (Yusef-Zadeh et al.||2007). Therefore, 
we conclude that LECKe scenario cannot be a major con¬ 
tributor to the remaining X-ray emission from Sgr B2 in 
2013. 


6.3. Cosmic Ray Ion Bombardment 

Although it is clear that the Sgr B2 Fe Ka emission 
decreases since 2001, The same level of the Fe Ka emis¬ 
sion in 2012 and 2013, and the lacking of measurements 
between 2005 and 2011 do not allow us to determine 
whether the Fe Ka emission was still decreasing in 2013 
or it has reached the constant background emission level. 
Future observa tions of Sgr B2 can distinguish between 
the two cases. Dogiel et al. (2009) estimated that CR 
protons could contribute to about 15% of the observed 
maximum Fe Ka flux obtained around 2000, although 
this prediction is highly model dependent. In 2012 and 
2013, the Sgr B2 Fe Ka emission reached 20% of that 
measured in 2001 (Section 4), therefore the 2012 and 
2013 Fe Ka emission could be mainly due to CR proton 
bombardment. 

The 2013 NuSTAR Sgr B2 spectrum can be fit well 
with the self-consistent LECRp model. The required 
metallicity of ZjZ^ = 2.5j'ji[ Q overlaps with the GC en¬ 
vironment metallicity of one solar to twice solar. Assum¬ 
ing all the X-ray emission from the central 90" Sgr B2 
in 2013 is due to CR proton bombardment, the re¬ 
quired total power of the CR proton of energies between 
Emin = 10 MeV and E^ax = 1 GeV in the cloud re¬ 
gion is dW/dt = (1.0 ± 0.3) X 10^® erg s“^. Taking the 
uncertainty of Emin into account, the required GR pro¬ 
ton energy ranges in dW/dt = (0.4 — 2.3) x 10^^ erg s~ ^ 
(Section 5.1.4). According to |Tatischeff et al.| (|2012|), 
there is an additional 40% power comes from a-particles 
with CalCp = 0.1, the final required total kinetic GR 
ion power is (0.6 — 3.2) x 10^® erg s“^, which is roughly 
10% of the steady state mechanical power supp lied by su¬ 
pernovae in the inner ^ 200 pc of the Galaxy [Tatische^ 
et al.| (|2012D. The required power injection is quite sig- 
nihcant and and probably difficult to fully be accounted 
with typical sources. 

The power deposited into the cloud is lower than the 
incident CR ion power, because CR ions with energies 
lower than Emin cannot penetrate the cloud and that 
with a path length of A = 5 x 10^“^ cm“^ those ions with 
energies higher than 180 MeV ca n escape from the clou d 
without depositing energy in it (Tatischeff et al. 2012). 
Therefore, the power deposited by CR ions into central 
Sgr B2 is Wd ^ 8 x 10^® erg s“^. With the central 90" of 
Sgr B2 mass of M ~ (0.5 — 2) x 10 ®Mq estimated based 
on the simplified Sgr B2 densit y profile (Section 1) an d 
a total mass of M = 6 x 10 ®Mq (Lis & Goldsmith^l990), 
the corresponding ionization rate can be estimated to 
be Cff ^ (6 — 10) X 10~^® s“^ using Equation (11) 


Tatischeff et al. (2012). This CR ionization for the 


dense materials in Sgr B2 is comparable to the GC CR 
ionization rate of ~ (1~3) x 10“^® s“^, which is found 
to be uniform thr oughout the GC on scales of 200 pc (e.g. 
Goto et al.||2011|). 


In the following we compare the application of LEGRp 
model on Sgr B2 and the Arches region and discuss 
possible LECR ion sources. The X-ray emission from 


the Arches region was interpreted as GR ion bombard- 
ment before a signif icant decrease in Fe Ka was revealed 
(Clavel et al. 2014). In the LECRp scenario, the CR 
spectral index is constrained by the broadb and NuS- 


TAR Arches region spectrum to s = I.65lg ®| (Krivonos 
et al. 2014), which agrees with that derived from the 


Sgr B2 specrtum (s = 1.9^®'®). The GC LECR ion 
source is suggested to be Galactic supern ovae or star 
accretion onto Sgr A* |Dogiel et al. 2013 and the ref¬ 
erences therein). The arches spec trum requires an ion- 
ization rate {(h 10“^® s“®, Tatischeff et al.|2012 ) 

higher than that of Sgr B2 by a factor of ~ 10 — 20, 
significantly exceeding the CR ionization r ate in the GC 
environment. Recently, Clavel et al. (2014) found a 30% 
flux drop from the Arches cluster in 2012, which indi¬ 
cates that a large fraction of the non-thermal emission 
of the Arches is likely due to the past activity of Sgr A*. 
This could explain the high ionization required by the 
Arches spectrum, assuming all the X-ray emission is due 
to LECRp. If the CR ionization rate of the Arches clus¬ 
ter is comparable to the GC CR ionization rate, LECRp 
bombardment should contributes to only a few percent 



contribution to the X-ray emission. While it is uncer¬ 
tain when the Arches cluster X-ray emission will reach 
the constant background level, it is likely that Sgr B2 has 
reached or will soon reach the background X-ray emission 
level. Therefore, the remaining emission from Sgr B2 will 
be a unique powerful tool to probe the CR population in 
the GC in the X-ray energy bands. The LECRp scenario 
will stay valid as long as no further decrease is observed 
in the coming years. 


6.4. The Nature of GO. 66—0.13 

GO.66—0.13 is an elliptical feature with a major ra¬ 
dius of ~ 5 pc and a minor radius of ~ 3 pc. The 
illuminating source photon index (F = 1.4 ± 0.5) is 
compatible with that of the central region of Sgr B2 
(F = 2.2±0.4). The intrinsic column density of Nnii) = 
3.0j'j® g X 10^® cm“^ is lower than the column density of 
the central part of Sgr B2. The primary source luminos¬ 
ity based on its 2012 peak Fe Ka emission is Ls = (1.3 ± 
0.5) X 10®®(i?/100 pc)2 (Zpe)”^ erg s“®, or L3-79 keV = 
(4.3 ± 1.6) X 10®®(i?/ lOOpc)^ (Zpe)”® erg s“^. It is con¬ 
sistent with the primary source luminosity required by 
the central region of Sgr B2. The Fe Ka line intensity 
of GO.66—0.13 exhibits fast decay at 6.4 keV. The lin¬ 
ear fitting shows that the GO.66—0.13 Fe Ka emission 
demonstrates an increasing trend prior to 2013 (Fig¬ 
ure 5). However, the increase could be much sharper 
than what is shown in Figure 5. We do not have good 
enough statistics to distinguish between a linear increas¬ 
ing from 2000 to 2012 and a sharp peak in 2012 on top 
of a baseline Fe Ka emission. In 2012, the peak Fe Ka 
emission from GO.66—0.13 was brighter than the Sgr B2 
core Sgr B2(M). The non-detection of its Fe Ka emis¬ 
sion in 2013 suggests a fast decrease in Fe Ka line flux. 
The fast flux decrease with a life time of t ~ 1 year 
roughly matches the light crossing time of the peak 
within GO.66—0.13 (r/c ~ 2 years). 














































14 


S. Zhang et al. 


Both Sgr B2 and G.66—0.13 are about 100 pc away 
from Sgr A*, while the GO.66—0.13 X-ray flux reached 
the peak 12-18 years after Sgr B2. If G.66—0.13 was 
illuminated by the same Sgr A* outburst that illuminated 
Sgr B2, it should be 16-23 pc behind the Sgr B2 center 
in the line-of sight assu ming Sgr B2 is 13 0 pc in front 
of the projection plane (Reid et al. 2009). Therefore, 
GO.66—0.13 could be a inolecular clump, with a higher 
local density than its surrounding environment, located 
in the Sgr B2 envelope which extends to about 22.5 pc. 

7. CONCLUSIONS 


We studied the X-ray emission of the central 90" of 
Sgr B2 and a newly discovered cloud feature GO.66—0.13 
with NuSTAR. We present their broadband spectra, 
morphology and time variability and discuss two possi¬ 
ble origins of the observed X-ray emission in 2013: X-ray 
reflection nebula and cosmic ray bombardment. 

• Central 90" of Sgr B2: The substructures of Sgr B2 
are resolved at sub-arcminute scales at hard X-ray en¬ 
ergy (> 10 keV). The hard X-ray continuum emission 
in 10-40 keV reveals two compact star-forming cores 
Sgr B2(M) and Sgr B2(N) surrounded by diffuse emis¬ 
sion with the western side brighter than the eastern side. 
The central 90" region is marginally optically thin on 
average. 

• Main Compact Cores Sgr B2(M) and Sgr B2(N): 
Gompact cores Sgr B2(M) and Sgr B2(N) are resolved 
above 10 keV for the first time. Both the 6.4 keV line 
emission and the 10-40 keV continuum emission peak at 
the location of the main compact core Sgr B2(M). While 
Sgr B2(N) is significant in 10-40 keV continuum emis¬ 
sion, it is not detected at 6.4 keV in 2013 perhaps due 
to higher local absorption. The Sgr B2(M) spectrum re¬ 
quires a high optical depth of tt = l.dl^;^, consistent 
with the Herschel measurement, suggesting it is not op¬ 
tically thin. 

• XRN vs. LECR: After a decade of decreasing, the 
2012-2013 Fe Ka emission reached 20% of that in 2001, 
but remained at the same level during 2012-2013. The 
lack of data between 2005 and 2011 does not allow us 
to determine whether the 2013 Fe Ka emission kept de¬ 
creasing (which could be explained by the XRN or the 
LEGRe model) or had reached the constant background 
level (for which the LEGRp model is favored). We first 
excluded the LEGRe scenario based on unphysical best- 
fit parameters, no matter whether the Fe Ka emission is 
decreasing or not. The significant 6.4 keV line of Sgr B2 
in 2013 requires the Fe abundance to be at least 3.4 solar, 
significantly exceeding the current measurements of the 
GG metallicity. The best-fit minimum electron energy 
is far below 100 keV. With such low energies, the non- 
thermal electrons are not able to penetrate Sgr B2 even if 
there is a local GR particle acceleration site. Therefore, 
we conclude that the LEGRe scenario cannot be a major 
contributor to the remaining level of the X-ray emission 
from Sgr B2 in 2013. 

The 2013 Sgr B2 X-ray emission can be best explained 
by the XRN scenario if the X-ray emission is still de¬ 
creasing. We examine the XRN scenario with an ad hoc 
XRN model and a self-consistent XRN model. Due to 
lack of geometrical information in this ad hoc model, in¬ 
terpretation of the intrinsic column density of A^h(*) = 
(5.0 ± 1.2) X 10^^ cm“^ calls for caution. With the self- 


consistent XRN model MYTorus, we derive the intrinsic 
equatorial column density to be Abj(i) = (1.0 ± 0.4) x 
10^^ cm“^, twice that derived from the ad hoc model. 
This corresponds to tt = 0.67±0.27, confirming that the 
central 90" region is marginally optically thin on average. 
In the XRN scenario, we develope an analytical method 
to calculate the primary source luminosity via Comp¬ 
ton scattering with the 10-40 keV scattered continuum 
(Equation 5). With this approach the primary source lu¬ 
minosity in any X-ray energy band can be derived with 
the observed continuum emission. The resulting pri¬ 
mary source luminosity is L 3-79 kev{scat) = (h.ot^'g) X 
10^®(d/100 pc)^ erg s“^. This is remarkably consistent 
with the primary source luminosity calculated via the 
photoelectric absorption process based on the updated 
source spectral shape, resulting in ^3-79 kev(cibs) = 
(5.0 ± 2.3) X 10^®(d/100 pc)^ erg s“^. We find that 
Ls{abs)/Ls{scat) is of order of unity, confirming the self- 
consistency of the XRN model. 

In case the Ee Ka emission has reached the back¬ 
ground level in 2013, the reflected X-rays could have 
completely faded and the LEGRp process could be a 
major contributor. The required total GR ion power is 
dW/dt = (0.6 - 3.2) X 1039 erg s"!, about 10% of the 
mechanical power supplied by supernovae in the inner 
^ 200 pc of the Galaxy. The CR ionization rate is found 
to be Ch (6 — 10) X 10“ 35 H“3 s“3, consistent with the 
GR ionization rate in the GG environment. If the Sgr B2 
X-ray emission has indeed reached the background level, 
it would be a powerful tool to constrain the CR ion pop¬ 
ulation in the GG. 

• The new cloud feature CO.66—0.13: The Fe Ka 
line flux of GO.66—0.13 reached its peak in 2012 and 
quickly diminished within ~ 1 year. The fast variabil¬ 
ity is best explained by the XRN scenario. The required 
primary source luminosity L3-79 keV = (4.3 ± 1.6) x 
103®(d/100pc)9Zpg erg s“3 is consistent with that de¬ 
rived from central Sgr B2. Assuming G.66—0.13 is illu¬ 
minated by the same Sgr A* outburst that illuminates 
Sgr B2, it should be 16-23 pc behind the Sgr B2 center 
along the line-of sight, which is still within the Sgr B2 
region. GO.66—0.13 could be a molecular clump with a 
local column density higher than surroundings located in 
the Sgr B2 envelope. 

• Other Calactic center molecular clouds and Sgr A* 
outburst history: In the NuSTAR Galactic plane sur¬ 
vey, other Galactic center giant molecular clouds were 
also detected well above 10 keV, including the “Bridge” 
and GO. 13—0.13. The same methodology to derive the 
primary source luminosity is applied to each cloud and 
reconstruct the Sgr A* outburst history (Mori et al. sub¬ 
mitted). 
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of the torus 


Figure 6. A spherical cloud (red circle) in a virtual torus (grey area) shows how the MyTorus model geometry is related to the GCMC 
geometry. The illuminating source (Sgr A* for GCMC) is shown by the yellow star. In the MYTorus model, the incident angle Oohs is defined 
as the angle between the symmetry axis of the torus (dotted grey line) and the LOS (black vertical line). The equatorial column density 
Nu is defined over the minor diameter of the torus. When applying to the GCMC case, the quasi-spherical cloud can be considered as part 
of the torus. Incoming X-ray photons from the illuminating source are scattered by the cloud at an angle 0 to the direction of the observer. 
The green lines show the travel paths of the X-ray photons, dproj is the projected distance between the cloud and the illuminating source 
as seen by the observer, and dios is the LOS distance between the cloud and the projected plane. A face-on view of the torus (0o6s = 0°) 
corresponds to the case where the cloud locates in the projected plane noted by a dashed black horizontal line {9 = 90°, dios = 0). The 
figure is from Appendix B1 in Mori et al. (submitted). 
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APPENDIX 

APPLICABILITY OF MYTORUS MODEL TO SGR B2 SPECTROSCOPY 

The MYTorus model was originally developed to study X-ray spectra of Compton-thick AGNs. The assumed geometry 
is a torus with a uniform distribution of neutral material reflecting incoming X-rays from an illuminating source at the 
center. Three model components are offered: the transmitted continuum (MYTZ), the scattered continuum (MYTS) 
and Fe fluorescence lines (MYTL). The model allows a range of values for three key model parameters: the illuminating 
source photon index P = 1.4 — 2.6, the equatorial hydrogen column density (corresponding to the minor diameter of 
the torus) Nh = 10^^ — 10^® cm“^ and the inclination angle dobs = 0° — 90°. Gonsidering a spherical molecular cloud 
as part of a virtual torus, the MYTorus model can be applied to molecular cloud spectra with some limitations. Figure 
A6 shows the geometry of a molecular cloud and a virtual torus. In the MYTorus model, the observation angle dobs 
is defined as the angle between the light-of-sight (LOS) and the symmetry axis of the torus. In contrast, most GC 
molecular cloud publications use the scattering angle d. A face-on case in the MYTorus model {dobs = 0°) corresponds 
to a cloud in the same projected plane as the illuminating source, where the scattering angle is 6* = 90°. For three 
key assumptions of the MYTorus model, we discuss the valid parameter space where the model is applicable to the 
molecular cloud X-ray reflection spectra in the following. 

Toroidal Reflector 

In the face-on case {dobs = 0°), MYTorus gives an accurate solution for a fully illuminated quasi-spherical molecular 
cloud, as different azimuthal parts of the torus all scatter at the same angle d = 90°. To derive the incoming X-ray flux 
that illuminates the cloud (red circle in Figure A6, which is part of the grey virtual torus), we simply need to rescale it 
by the solid angle ratio of the torus (fixed to Q/2n = 0.5) and the cloud. For Sgr B2, its solid angle respect to Sgr A* 
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is ~ 2.89 X 10“^. When the inclination angle deviates from 6obs = 0°, i.e. the face-on case, different azimuthal parts 
of the torus scatter incoming X-rays at different angles. The reflected spectrum thus shows variation and becomes 
inaccurate for a quasi-spherical molecular cloud. However, the scattered component MYTS does not vary strongly 
with 9obs as long as 9obs ^ 60° and Nh < 10^^ cm“^ (See Figure 17 in Mori et al. 2015). When 9obs > 60°, the 
spectral variation starts to become significant, for the following two reasons. Firstly, as the inclination angle increases, 
the X-ray photons back-scattered by the side of the torus hit the closer side of the torus before reaching the observer, 
thus are subject to further absorption. Secondly, multi-scattering can become important and cause angular-dependent 
X-ray flux, although it is negligible at Nh ^ 10^“^ cm“^. 

Next we determine the systematic error of measuring Nh, F and model normalization based on the MYTorus model. 
We simulated MYTS spectra for 9 = 60° and Nh = 10^^ — 10^'^ cm“^ and fit with the MYTS model with 9obs set to 
0°. The deviations of the best-fit value of these model parameters from their input value are adopted as the systematic 
errors. At Nh = 10^^ cm“^, the deviation of Nh, F and model normalization from their input are 10%, 1% and 7%, 
while at Nh = 10^^ cm“^, the deviations increase to 25%, 3% and 10%. Therefore we conclude that the reflected X-ray 
spectrum model MYTorus is not sensitive to the reflector geometry in the NuSTAR energy band as long as 9obs 60° 
and Nh ^ 10^"* cm“^. Similar to the face-on case, the incident X-ray flux can be derived by rescaling it with the 
solid angle ratio of the torus H/Ttf = 0.5 and central 90" of Sgr B2 cloud H/Ttf = 2.89 x lO”'^ x (100 pc/R)^ with 
systematic error < 10%. However, the conversion from the incident angle in the MYTorus model to the position of the 
cloud cannot be well established. We note that the Nh measured by the MYTorus model is an averaged value over 
the torus: it does not i nto account the specific geometry of the studied cloud and its possible partial illumination (see 
e.g. Odaka et al.||2011 ). The systematic errors are estimated solely based on the MYTorus model and might therefore 
be underestimated, 'ine more sophisticated model under development will better constrain the Nh (Walls et al. in 
prep). 


Fe Abundance Fixed to Solar 

The Fe abundance is important in determining the incoming X-ray flux solely from measurements of Fe fluorescent 
lines, while the 3-79 keV broadband spectrum of NuSTAR also provides the energy range where Compton scattering 
dominates (> 10 keV) over Fe fluorescence. The incoming X-ray flux is determined consistently from both the Compton 
scattering process with the MYTS component and the Fe fluorescence process via the MYTL component. The best-fit 
model parameters {Nh, F and normalization) of MYTS do not vary significantly when the MYTS model is fit to the 
cloud X-ray spectra with or without the 6-10 keV energy range which covers the 6.4 keV Fe Ka line, the 7.06 keV Fe 
K/3 line and the 7.1 keV Fe K edge. Data tables for different Fe abundances in the range of Zpe = 0.5 — 3.0 will be 
implemented to the modified XRN model for molecular clouds. 

Uniform Density 

The incoming X-ray flux from the illuminating source will not be significantly affected by non-uniform density profiles 
as long as Nh < 10^"^ cm“^, where multi-scattering is negligible. If the density is not uniform, large clouds can include 
very dense clumps where multi-scattering effects can be significant. For Sgr B2, the largest optical depth is ~ 1 
in the direction of the densest cores Sgr B2(M), in which case the multi-scattering effects are negligible. The current 
XRN model does not address the complicated density profile of Sgr B2, which contains compact cores, clumps and 
an overall decreasing density profile. A more sophisticated XRN model with a reasonable density profile of Sgr B2 
implemented is under development (Walls et al. in prep) and will reduce the uncertainties caused by non-uniformity 
of Sgr B2. 
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